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Abstract: We compute the first perturbative correction to the static correlation lengths 
of light quark bilinears in hot QCD with finite quark chemical potentials. The correction is 
small and positive, with ^-dependence depending on the relative sign of chemical potentials 
and the number of dynamical flavors. The computation is carried out using a three- 
dimensional effective theory for the lowest fermionic Matsubara mode. We also compute 
the full correlator in free theory and find a rather complicated general //-dependence at 
shorter distances. Finally, rough comparisons with lattice simulations are discussed. 
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1. Introduction 

With the ongoing and future heavy ion experiments at RHIC and LHC, we have at hand 
a large body of experimental data on QCD at extreme conditions. However, comparing 
this data with the theory predictions is far from trivial, because we cannot usually reliably 
compute the measured physical quantities directly from QCD, even in the case when the 
plasma reaches an equilibrium. 

The only systematic non-perturbative approach is lattice gauge field theory in four di- 
mensions. With an increasing amount of computer power available, it has become possible 
to determine many quantities starting from first principles, while avoiding the problems 
arising from the strong coupling and higher order infrared divergences that make per- 
turbation theory so hard near the phase transition. Lattice simulations, however, have 
potentially large systematic errors arising from light dynamical quark degrees of freedom. 
In the presence of chemical potential, which is the situation studied in this paper, the 
situation is even worse, since the quark determinant becomes complex and conventional 
importance sampling does not work. 
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There is clearly room for other methods, whose results, even if they have their own 
limitations, could be compared with the lattice to gain insight on the shortcomings of both. 
Perturbation theory suffers from infrared problems at high orders [1], but these are only 
related to static, soft gluons. This has led to a technique called dimensional reduction, 
where the nonperturbative infrared behavior is pushed to a simpler effective theory for soft 
modes only, allowing the lower order contributions to be computed reliably in perturbation 
theory. Of course, at moderate temperatures the QCD coupling is large and higher order 
terms cannot be neglected, but in this way computationally intensive nonperturbative 
methods can be applied to a simpler 3d theory, while the terms with dynamical quarks 
and quark chemical potentials can be computed perturbatively in four dimensions. On the 
other hand, at very high temperatures the perturbation theory result computed in this way 
should be reliable because of the asymptotic freedom. 

In the bosonic sector dimensional reduction has been used successfully in both QCD 
and electroweak theory to compute quantities like pressure [2-8] and various gluonic cor- 
relators [9-14]. Fermionic modes, on the other hand, are usually integrated out from the 
effective theory and only affect the parameters of the dimensionally reduced theory. How- 
ever, there are many interesting observables constructed of quark fields which are sensitive 
to infrared physics and whose perturbative treatment thus requires some kind of resumma- 
tions, conveniently organized as successive effective theories. In this paper we concentrate 
on the long distance behavior of correlators of mesonic operators composed of light (mass- 
less) quark flavors [15]. 

The first efforts on treating this kind of correlators in a dimensionally reduced theory 
date back almost 15 years [16,17]. However, those works did not systematically include all 
terms of a given order. In [18] it was shown that the dimensionally reduced theory can be 
formulated in terms of nonrelativistic quarks, and the correct power counting of different 
operators was established. With similar methods we [19] were able to consistently compute 
the 0(g 2 ) corrections to the screening masses of mesonic operators. The current paper aims 
to extend that work for situations with nonzero quark chemical potentials. These masses 
have been measured recently on lattice both at zero [20,21] and finite density [22]. Near 
the QCD phase transition these correlators have also been studied at finite temperature 
and density using the Nambu-Jona-Lasinio model [23] . 

This paper is organized as follows. In Sec. |2|we define the correlators that we intend to 
study and compute the leading order (free theory) results. The derivation of the effective 
theory along the lines of [19] is carried out in Sec. ||, as well as the matching of the pa- 
rameters to the level required for computing the screening masses. In Sec. Q we formulate 
the problem in the context of the effective theory and solve the screening mass numeri- 
cally. Discussion of the results and comparison with lattice data can be found in their 
respective subsections. Our conclusions are presented in Sec. || Appendices [A| and || con- 
tain the details of computing the leading order correlators and quark-antiquark potential, 
respectively. 
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2. Static correlators in high temperature QCD 



The quark sector in finite temperature QCD is described by the Euclidean Lagrangian 

Cl = ^D fl + M)iP, (2.1) 

where the covariant derivative of quark field is defined as 

= - igA^T a 4>, (2.2) 

with T a the generators of the fundamental representation of SU(iV c ). The quark field ip 
is an iVi?-component vector in flavor space, and the Np x Np mass matrix M is taken to 
vanish in most of what follows. 

Finite quark densities can be introduced via chemical potentials /if. Inserting /if into 
the path integral leads to the effective Lagrangian (with M = 0) 

A = E ^(^^f + W7o)Vr • (2-3) 
f 

In perturbation theory the effect of the chemical potential is then to give a constant imag- 
inary part to the momentum zero-component, po — > po — ifif. The chemical potentials for 
each quark flavor can be chosen independently when no weak interactions are present, while 
weak interactions would only conserve certain specific combinations of baryon and lepton 
numbers, see e.g. [24]. We try to give our results for general /if when possible, but will 
mainly use isoscalar (/x u = /id = A*s) an d isovector (/x u = — /id = fJ-v) chemical potentials 
to discuss the effects of various combinations of /if. 

With the quark fields we can define various bilinear objects of different spin and flavor 
structure, 

O a = ^F a T^, (2.4) 

where V is one of {1,75,7^,7^75} for scalar, pseudoscalar, vector and axial vector objects 
S a , P a , V® and A^, respectively. The identity matrix F s and the traceless matrices F a 
generate the flavor basis, 

F s = l Np , Ti[F a F b ] = ^S ab , a,b=l,...,N F -l. (2.5) 
We are interested in correlators of the form 

C q [O a , O b ) = J^dr J d 3 x e iq x (O a (r, x)O 6 (0, 0)), (2.6) 

or, in configuration space, 

C*[O a ,O b ]^ / 1/T dr(O a (r,x)O b (0,0)). (2.7) 
Jo 

Because of rotational invariance, the general structure of correlators in terms of the 
three-momentum q and the chemical potentials /if, fi g of the quark flavors in the operator is 

C n [S a ,S b ], C n [V a ,V b ] = J2 F ?s F i /(^W^g), (2-8) 

fg 

W>*?1 = E^gf [(*> " q ~f) % 2 ^f,%) + q -fKq 2 ,M^ g ) 
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(a) (b) (c) 

Figure 1: Diagrams contributing to the meson correlator. 



and similarly for P a and A^. Here q = |q|. We expect that the functions /, t and / are 
dominated by simple poles near the origin q 2 = 0, corresponding to a set of bound states 
in a 2 + 1-dimensional theory obtained by analytic continuation. Their Fourier transforms 
then exhibit exponential falloff at large distances, and the coefficients of this falloff are 
referred to as screening masses. 

Without a loss of generality we can select the direction of x as the X3-direction. We 
make this choice and further average over the xiX2-plane, giving a more easily calculable 
correlator 



C z [0\O b ] = f d 2 x ± <7(xx,*)[0 o ,0 6 ] = j'^dT f d 2 x ± (O>,x ± ,z)O fe (0,0,0)>. ( 2 



10) 



2.1 Leading order in perturbation theory 

At very high temperatures the strong coupling constant approaches zero because of asymp- 
totic freedom. Perturbation theory should then be applicable, and we can expand the result 
in powers of a s around the perturbative vacuum. The leading order result for the correlator 
in Eq. ( 2.10| ) is given by the free theory one loop diagram of Fig. |](a). 



For isoscalar chemical potential the leading order correlator can be computed exactly. 
A straightforward calculation, the details of which can be found in Appendix [5j gives 

NT/ 1 \ 

C z [S a ,S b ] = 5 ab \ 2ttT coth 27rTz cos + 2u sin 2pz + -cos2//z] (2.11) 

8irz sinh 2-kTz V z 



8vr 2 z 3 

xab ^cT 2 1 -2wTz 

2 z 



1 \ u 
1 + J cos 2yz + — sin 2/j,z 



+ 0(e 



-AttTz\ 



which for [i = agrees with the previous zero density computation by Florkowski and 
Friman [25]. In Fig. ^ we have plotted the expression in Eq. ( |2.11[ ) for various \i. The 
wavelength of the oscillation from terms like cos(2/iz) is L = n/fi, while the exponential 
falloff is characterized by the screening length = 2vrT. When / M > or y. < 2tt 2 T, the 
oscillations become important at distances where the correlator is exponentially suppressed 
by screening and show only as small ripples in the tail part of the C z . For larger chemical 
potentials there are strong oscillations inside the exponentially decaying envelope. 

This behavior is verified by Fig. ^, where the oscillatory behavior cannot really be 
seen for fj, < 0.57rT, but the correlator just falls off faster than for zero density. At short 
distance this can be viewed as increased screening, but we will define the screening mass 
strictly as the coefficient of the exponential falloff, or the real part of the pole location. 
Nevertheless, the imaginary parts will be computed as well. 
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The free correlator with arbitrary chemical potentials is much harder to compute. In 
Appendix [A| we compute an approximate correlator by restricting the sum over Matsubara 
modes to u>n = ±7rT, since they dominate the long distance behavior of the correlator. 
The result can be given in terms of the exponential integral function Ei(x), and to leading 
order in 1/z it is 



,N r T 



C z [S a ,S b ] = Y^F^Ff^- (l - e - A ^ T ) -e-^ Tz 



1 



1 



A/iy 



A Wj (2vrT) 2 + p* 



A Wj 



A/iy (2vrT) 2 + p\ 



2nT cos /iy z 



Mij sm^ijZ 



(2.12) 



where A/Uy = pb\ — /ij and /iy = /i; + //j. It should be noted that at this order the scale of 
the oscillation is p^, which vanishes for isovector chemical potential. The subleading terms 



contain oscillations with A/zy, as can be seen in Eq. (A. 3). 

Pseudoscalar correlators differ only by a sign from Eqs. ( fPH) and The (axial) 

vector correlators can be likewise computed in the free theory and have the same general 
behavior C z ~ cos(pijz) exp(— 2ttTz), while the structure is a bit more complicated. 



3. Effective theory 

We would like to compute the next-to-leading order correction to the screening mass in 
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perturbation theory. As usual in finite temperature field theory, simply evaluating the 
next-to-leading order diagrams in Fig. |](b,c) is n °t enough, but we have instead an infinite 
set of diagrams all contributing at the same order. 

First of all, the temporal gluonic zero Matsubara mode ^4q ^ gets a self-energy correction 
of order g 2 T 2 , which for small momenta should be included in the propagator as an electric 
mass. The introduction of thermal mass resums the diagrams with an arbitrary number 
of self-energy insertions on the gluon line, and it is required in order to get rid of the 
infrared divergences that the low momentum electrical gluon modes give rise to in naive 
perturbation theory. 

Since we are computing corrections to a mesonic operator near the threshold of pro- 
ducing two free quarks, the quarks can be almost on-shell, with ~ 0(l/g 2 T), which 
compensates for the factors of coupling from the vertices. The diagrams with an arbitrary 
number of low momentum gluonic zero modes exchanged between quarks should therefore 
be resummed to get a consistent first order correction to the propagation of free quarks. On 
the other hand, we do not need to compute the gluon exchange diagrams like in Fig. |(b) 
for non-zero gluon modes, since only the one-gluon diagram would contribute at order g 2 , 
and this diagram alone cannot change the pole location, i. e. the screening mass. Of course, 
if we were computing the overall normalization or the functional form of the correlator at 
short distances, then these diagrams would contribute. 

A convenient way of organizing these resummations is the effective field theory method. 
The goal is to write a simpler theory for those modes that require special treatment, namely 
the gluon zero modes and the lowest fermionic modes. These modes reproduce the correct 
infrared behavior of gauge theory at high temperature [26], the relevant scales being of the 
order gT and g 2 T, while the contribution of all the other modes can be included in the 
parameters of the theory. With the sum over Matsubara modes suppressed, the remaining 
theory is three-dimensional. 

In the bosonic sector the dimensional reduction to the physics of gluonic zero modes 
is well-known. At finite density the Lagrangian has the form [13] 

£ c b ff = ±TrF 2 + Tr[A, A ] 2 + m 2 E TrA 2 + ^ E ^ A o + X e {^ A lf + ^< (3-1) 

f 

with i = 1,2,3, Fij = ig^ 1 [Di, Dj], Di = di — ig^Ai and the parameters to the order we 
need them are 



ml 



^ 2 (f + ^ + 9 2 E =9 2 T, \^=0{g±T). (3.2) 



The cubic and quartic Aq interactions can be ignored, since they contribute to meson 



correlators only at order g° and higher. The theory in Eq. (3T) describes a 3-dimensional 
gauge theory with a massive adjoint scalar Aq. 

The fermionic sector consists of the lowest fermionic modes u = ±7rT. The reason we 
are leaving out the other modes is not that they would be substantially heavier (irT <C 3irT 
is not really true), but because the lowest modes dominate the correlator at large distances. 
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We are expanding around a state consisting of a free quark-antiquark pair, and the relevant 
expansion parameter is the "off-shellness" p$ + p\ + p|, which near the screening pole 
y>3 ~ ±i(irT — i/i) is of the order of the dynamical scales gT,g 2 T. 

In the three-dimensional theory the lowest fermionic modes can be viewed as massive 
particles of mass \uj±i\ = irT. As noted above, the typical momenta are much smaller than 
this. The correlators then appear as bound states of massive quarks in (2+l)-dimensional 
theory, a situation resemblant of quarkonia in four dimensions. At T = the properties 
of quarkonia can be studied with the help of an effective theory known as non-relativistic 
QCD (NRQCD) [27], and the fermionic sector of our final effective 3d theory Eq. ( 3. 12| ) 



looks very similar to this theory, although the physical interpretation of the parameters is 
quite different. Because of similar appearance, we label our effective theory NRQCD 3 . Its 
derivation will be presented below. 

3.1 Tree-level NRQCD 3 

The dimensional reduction in the fermionic sector using NRQCD methods was first studied 
by Huang and Lissia in [18], where they also computed a number of one-loop corrections 
to the parameters. A similar derivation in a slightly different formalism was presented in 
our previous work [19]. We will here review the reduction in the case of finite chemical 
potentials, and point out the differences. 

For a field VK X ) corresponding to the quark Matsubara mode u n the Euclidean La- 
grangian (|2.3| ) reads 

£ q = ip [z7 w„ + 70^ - ig"foAo + 7 fe -D fc + 73^3] i>, (3.3) 

where k = 1,2 and Aq is the gluonic zero mode. Note that the interaction vertex with Aq 
does not mix the different fermionic modes, so we have separate terms like Eq. fl3.3j ) for 
each mode. We have used the rotational invariance to measure the correlations always in 
the a^-direction, and therefore separate the term from the transverse ones. 

In the nonrelativistic region we should be able to treat the particle and the antiparticle 
separately. We choose a nonstandard representation for the Euclidean Dirac matrices, 

70H—1, *-(t4J' *-•«<•'>• <3 - 4) 





so that 



eijPiVj (Po ~ iP3)h ) ' 



(3.5) 

where is antisymmetric and 612 = +1- This is diagonal in pq and P3, so writing 
gives us 

£q = iX ] (Po ~ gA + As) X + (Po ~ gA - D 3 )(p + ^e^Diajx - x^ijDi^j^, (3-7) 
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where po again includes the chemical potential, po = iv n — ifj,. 

As stated above, we are expanding around the screening poles p 3 = ±ipo, where one 
component is light while the other one is heavy compared with the dynamical scales. Work- 
ing at the tree level, we can then solve the equations of motion for the heavy component 
and expand in 1/po, getting 



/>[> - fl- 1(1 + A; - — (D± + ^[°:-"M\, 



X 



po - gA - D 3 - — (p\ + ^[a h o^Fy 



+ 0(1). (3.8) 



At classical level the effective Lagrangian for quark Matsubara modes then consists of 
two nonrelativistic particles coupled to gluons and to a massive adjoint scalar, all in three 
dimensions. This can be rotated to a (2+l)-dimensional Minkowski space as in [18] to have 
a form similar to NRQCD, but for the screening states this is an unnecessary complication, 
and we choose to stay in Euclidean space instead. 

3.2 Matching at g 2 level 

When going beyond the tree level, the parameters in Eq. ( |3.8[ ) are subject to quantum 
corrections. We also need to consider all other possible operators allowed by symmetries 
but not present at the classical level. To keep the number of possibilities finite, it is 
necessary to set up a power counting for the possible operators. 

By definition, the effective theory contains quarks with |pj_| < gT and the off-shellness 
Ap3 = P3 ± ipo of the same order. Requiring the action to be of order unity, we have 

J dx 3 d 2 x ±X %X ~ 1, => X ~ V|x±| ~ gT (3.9) 
J dx 3 d 2 x ± AdlA ~ 1 A ~ (x 3 /x 2 ± ) 1/2 ~ g 1 / 2 ^ 2 . (3.10) 

For relativistic on-shell gluons the energy and momentum are of the same order, p 3 ~ p^. 
The quarks, on the other hand, are nonrelativistic, so their kinetic energy (p 3 in the 2+1- 
dimensional theory) is proportional to the momentum squared. As can be verified from 
the poles in quark propagators Eqs. ( 3.18| ),( |3,19| ), for nearly on-shell quarks in NRQCD 3 , 



with transverse momentum < gT, the off-shellness is 

Ap 3 ~ p 2 ±/po ~ g 2 T d 3 ~ g 2 T acting on quarks. (3-H) 

These estimates enable us to limit the number of possible operators. For example, 
any four-quark operator would be of the order g 4 T 4 , while the leading order behaves as 
g 2 T s . This means that in order to produce at most 0(g 2 ) corrections the coefficient 
would have to be of order g°, but this term is not present at tree level and can therefore 
be excluded. Likewise, there is a known kinematic correction D^/SpQ coming from the 
correct normalization of the spinors, but this is also of higher order (g 4 ) than we need in 
our calculation. 
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P p — q 



Figure 3: One-loop quark self-energy correction. 



The same power counting also shows that g E A ~ g 3 ^ 2 T, which is of higher order than 
the derivative term di in the transversal covariant derivative. Since d\ is already of order 
g 2 , we can leave out the transverse gluons. Similarly, the only term that needs matching 
beyond tree-level is the zero-point energy pq, which we will denote by M from now on. To 
order g 2 it is then sufficient to use the fermionic Lagrangian 

4s = iX ] (m - g E A + £> 3 - |T) x + itf (m - g E A - D 3 - 4>. (3.12) 

The one-loop correction to the parameter M can be computed by matching Green's 
functions computed in QCD and NRQCD 3 . To avoid the factors of g 2 arising from the 
normalization of the fields, we choose to match the location of the pole in quark propagator 
in the p3-plane. In both the effective theory and the full QCD this is a gauge invariant 
quantity, which corresponds to the quark pole mass in the (2 + l)-dimensional theory. 

We will use the gluon propagator with an infrared regulator A as well as a gauge param- 
eter £ and use dimensional regularization throughout to control the ultraviolet behavior. 
We will see that in the final results the dependence on £ vanishes and that we are allowed 
to take the limit A — ► 0, e — > 0. The gluon propagator in full QCD reads 



{A a Jp)A b M)=8 ab {2-K)^5^{p + q) 



p2 J p2 _|_ ^2 p2 p2_|_£^2 



(3.13) 



The one-loop self-energy correction to the quark propagator is given by Fig. ||, 

lq ,p- q ?( q +x 2 )- i9C *y q 

where we have used the convenient short-hand notation 



„ , I" J J (2tt) 



for dimensionally regularized sum-integrals in the MS scheme. 

The combination ip — S(p) should annihilate the on-shell spinor u(p) order by order in 
perturbation theory. When computing the order g 2 correction we can use the leading 
order results p 2 = and pu(p) = to drop terms proportional to p 2 and p, since £(p) is 
already multiplied by g 2 . The longitudinal part then vanishes and the result is independent 
of the gauge parameter £. The remaining expression is still rotationally invariant, but we 
can choose to match the Green's functions at some specific momentum and choose pj_ = 0. 
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When multiplied from the left by 70 , the expression becomes diagonal in the representation 
of Dirac matrices given by Eq. ( [3.4D , and we can concentrate on, say, pyoEjn, which is 



h^M = -x{D-2)g C F ^ {p _ q)2{q2 + x2y (3.16) 

We again use the tree-level relation p% = ip$ inside the integral, which then reduces to a 
sum of a bosonic and a fermionic tadpole, both infrared finite, giving the pole position 



P3 



(3.17) 



The NRQCD3 Feynman rules can be seen from Lagrangian in Eq. ( |3,12| ). In particular, 
the propagators are 

(Xu(p)xUq)) = Suv(2tt) 3 5(p - q) ~] (3.18) 

M + ip 3 + pj_/2p 

(MtiW) = S uv (27rf 5(p - q) , ~ l (3.19) 

M -ip 3 + pj_/2p 

or, in the configuration space, 

(xu(x)xUy)) = -id uv 9Mx3 - y3)) TT ^ — -«fe- J /3)-^-T 3 ; ) (3 20) 

2vr(a;3 - y 3 ) 

(MzWviv)) = -^K(y 3 - — e - M fc-^-t«!'. (3.21) 

27r(y 3 - x 3 ) 

These equations show clearly that for uj n > the x field propagates forward and <fi backward 
in X3, the time coordinate in the (2+1) -dimensional theory. For u n < the roles are 
reversed. 

In [19] we expanded the NRQCD 3 propagators in l/po> treating the kinetic term as a 
perturbation. This works well in the case of a single heavy quark, but, as discussed in [28], 
will run into problems in NRQCD at two-loop level. We will in this paper follow [28] in 
keeping the p^_/2po term resummed into the propagator. This also means that we have 
to use the multipole expansion, which in practice means that the gluons are not allowed 
to transfer transverse momentum to quarks. As a consistency check, we have shown that 
this formulation of the effective theory is able to reproduce exactly the free theory result 
Eq. (|2TT1) . 

On the effective theory side using the different propagators and adding chemical po- 
tentials does not change the result of [19]; the one-loop correction to quark self-energy 
vanishes at the pole in NRQCD 3 . The pole is then located at p 3 ~ iM and we can match 
these values to give 

M = Po + g 'C F ^- (l + -g,) _ u. - i» + fO^^T) (l + -g,) (3.22) 

The above equation was derived around the pole p% = ipo and for the (ll)-component of 
S(p), so it gives the mass of the x field. A similar computation for the pole p 3 = —ipo 
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shows that the <j> field has the same mass. We made no assumptions about the sign of u n 
here, and hence Eq. ( p. 22 ) is valid for all Matsubara modes. However, as we will see below, 



for negative modes it is — M that enters the equations as a mass paramater. It should be 
noted that for the lightest modes tu n = ±7rT, which eventually are the only important ones 
here, the real part of the mass is independent of fj,. 

4. Solution for the screening states 

With the effective Lagrangian of Eqs. (|]|), (fU|) at hand, we are now ready to compute 
the correlators. The various bilinears in Eq. ( |2.4| ) can be written as products of the x and 
4> fields with a 1 and flavor matrices between them. The details can be found in [19]. Since 
the propagators ( |3.18| ), ( |3,19| ) as well as the couplings in Eq. ( |3.12| ) are diagonal in spin 
indices, and (a 1 ) 2 = 1, the spin structure does not affect the result. The flavor structure, 
on the other hand, will be significant if we allow for different chemical potential for each 
quark flavor. 

As discussed in the beginning of Sec. ^, and verified by the power counting rules above, 
we have to take into account all diagrams of type Fig. [l|(b,c) with an arbitrary number 
of soft gluons. The matching computation already took care of the hard (u) b n ^ 0) part 
of Fig. |l](c), and the non-static modes in Fig. [j](b) do not contribute to the pole location. 
In the nonrelativistic theory we can take the gluon exchange to be instantaneous, and 
compute the static potential for the cj)*x pair by integrating out the gauge fields. The 
screening states can then be solved from the familiar Schrodinger equation. 

It should be noted here that it is enough to compute the potential to order g\. To see 
this, notice that the general structure of the potential can be written as an expansion in 

V(r) ~gl\nr + gir + 0(gir 2 ). (4.1) 
The Schrodinger equation gives to lowest order 

L&„ V (r)~!&hir, (4.2) 
Po or z 



so 1/r ~ yPo#E ~ 9^i or 9E r ~ 9- The leading term g^lnr coming from the one-gluon 
exchange is then sufficient for finding the g 2 correction to the screening mass. It should be 
noted that this term already gives a confining potential, so there is no qualitative difference 
in leaving out the linear term. There might still be infrared divergences when computing 
the scale inside the logarithm, but as we will see, this term is completely finite in NRQCD 3 
and no further resummations are needed. 

The potential can be found by introducing a point-splitting in the correlator to put 
the quarks a finite distance r apart, then finding the Schrodinger equation satisfied by the 
correlator and in the end letting r — > 0. This calculation for [i = was carried out in [19], 
and adding the chemical potentials and using propagators with the transverse kinetic terms 
included does not change the result, which is 

V(r) = (in + 1E - K (m E r)) . (4.3) 
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Since the computation of the potential was not shown explicitly in [19], we give the details 
in App. ||. 

The correlator C(r,z) for a 0*Xj pair, where i,j are flavor indices, then satisfies the 
equation 

' _ " 1 C(r,z) =2N c 5(z)5\r), (4.4) 



-V r 2 + F(r) 



where poij is the reduced mass of the system, 

1 _ 1 



+ 



1 



(4.5) 



and the + and — signs stand for positive and negative Matsubara modes us n = ±ttT, 
respectively. The exponential falloff is determined by the eigenvalues of the Schrodinger 
equation 



±(M- l + M l 



1 



-V 2 r + V(r) 



= "ifuii^o, 



(4.6) 



±2poy 

with some ground state wave function. It should be noted that the parameters of the 
modes ±ttT are related by M_ = — M* and po- = ~Pq + , and that the potential is real, so 
the complex conjugate of Eq. Q4.6[) gives 



mm,- 



rut 



%!!,+• ( 4 - 7 ) 

Assuming the ground state is unique (up to normalization) , the full correlator of Eq. (|2.1C| ) 
which is the sum over both positive and negative modes, behaves as 

C z [0\O b ] « 

« J d 2 R (<P* + (R, z)F a X+ (R, z) X \ (O)F fe + (0)) + (x- (R, z)F a cf>- (R, z)<j>*_ (0)F b X - (0)> 
« Y, F ?i F ^ [ e " ij exp(-m fu u )ij z) + e~ ia v exp(-m^ Mj z)] 



^FyFj- 2cos[Im(m fu ii ii j)z - ay] exp[-Re(m fulljij )z], 



(4.8) 



where is an unimportant overall phase of the <f>*Xj correlator. The correlator is evidently 
real as it should be and behaves as exp(— mz) cos(fxz), compatible with the leading order 
term in Eq. ( p. 12] ). In particular, the screening mass is given by the real part Re(mf u n), 
while the imaginary part contributes to the oscillation with z. Because of the relation 
Eq. ( |4.7j ), we can restrict the following analysis to the fields corresponding to the positive 
Matsubara mode +ttT. 

The eigenvalues of the Schrodinger equation ( |4.6| ) have to be found numerically. We 
change to dimensionless variables 

g|C F 2p ij 



2tt 



P 



giving for a spherically symmetric state 



2tt m| 



*n = 0. 



(4.9) 



(4.10) 
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The allowed values for Eq are found by integrating Eq. (|4.10| ) out from small r and 
requiring square integrability. This is somewhat more complicated than in the p = 



case, since now both p and Eq are complex, so Eq. ( 4.10 ) actually represents two coupled 
equations for Re(^o) an d Im(^o)- Modified versions of the integration routines we used 
in [19] can be applied, but the values of Eq now have to be searched in the complex plane 
instead of the real axis. In addition, the result depends on the chemical potentials in a 
nontrivial way, and the numerical computation has to be carried out separately for each 
/x/ttT. 

4.1 Results and discussion 

For isoscalar chemical potential (/i u = p^ = Ps) the result of solving Eq. ( f4,10| ) numerically 
is plotted in Fig. |||. Shown is the location of the pole in the complex So-plane for < 
Ps/ttT < 1. For negative values of ps Eq. Q4.10 ) shows that the solution is Eq(—/j,) = 



Eq(h). As can be seen from the figure, the pole moves on a quadratic curve off the real 
axis and towards larger Re(_Eo)> except for the Np = result. In the Np = 3 case 
we also study the situation where only two of the quarks have finite chemical potentials, 
whereas the density of the third quark vanishes. This should correspond to the situation 
in the heavy ion collisions, where the two incoming nuclei have no net strangeness, but the 
strange quark mass is small compared to the temperature. For isovector chemical potential 
the parameter p above is real, so the eigenvalues stay on the real axis. 
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Figure 5: The real part of Eq, for isoscalar (left) and isovector (right) chemical potential. 



The screening mass depends on the real part of Eq through 

Re(mfiji) = Re(Mi + Mj) + g^R e (E ) = 2itT + g 2 T^- Q + Re(£ )) , (4.11) 

where only Eq depends on the chemical potentials. In Fig. || we plot Re(-Eb) vs. fi for 
both isoscalar and isovector chemical potentials. In the isoscalar case the screening mass 
increases slowly with fig for dynamical fermions, while there is a decrease in the Np = 
(quenched quarks) case. The derivative of the real part with respect to \i vanishes at fj, = 
for all chemical potentials, which follows directly from the symmetry Eq{—^j) = E^fi). 

In the case of isovector chemical potential the mass decreases as a function of [iy when 
Nf < 2. At Nf = 3 the /iy-dependence in |>o oc (l+/i 2 ) and m\ oc (6+Nf+3Nf£i 2 ) cancels 
exactly in p, leading to a /iy-independent mass, but this is a plain numerical coincidence 
with no deeper physical significance. With even more dynamical flavors the mass would 
increase with (iy, but we have not done any numerics in those cases. 

The numerical value of the correction is small compared to the leading order result. 



For Nf = 2, Eq. ( |4.1l|) gives numbers in the interval 
Re(m fuU ) « 2vrT + g\ x { 



' 0.227, ns/nT = 1.0 
0.187, /i/vrT = 0.0 (4.12) 
0.166, pv/ttT = 1.0 



At T ~ 2T C the effective coupling is estimated [29] to be g\/T ss 2.2, so the first 
correction is about 6-8%. However, because the gauge coupling is large, this does not 
guarantee that the perturbative series would converge rapidly. 

To sum up, the general behavior of the screening mass in the presence of dynamical 
fermions is that it increases with isoscalar chemical potential and slightly decreases with 
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the isovector one. Note, however, that the behavior of the correlator differs from a simple 
exponential, in particular at short distances. When the correlator is measured only at 
short distances, it falls off faster because of the cosine term (see Fig. |2|), showing an 
apparent increase in screening. To see this more specifically, we note that simply fitting 
an exponential of the form C exp(— mz), where C and m are the fitting parameters, to the 
plane-plane correlator oc exp(— Mz) cos(pz) gives 

(m-M)(m + M) 2 = p 2 (3m + M) => m « M ( 1 + - -O. (4.13) 



M 2 4M 4 , 

One should therefore be careful when extracting the screening mass from the correlator, 
since fitting a wrong kind of a function will falsely interpret the oscillation length as part of 
the exponential falloff. This may also partly explain the considerably weaker p dependence 
observed on lattice in the isovector case. As can be seen in Eq. ( EU2l) , there is no oscillation 
at leading order when p u = — /id- 

It is not very clear how the screening mass should be determined in the finite density 
case if the functional form of the correlator is not known. At p = one often defines the 
effective (distance-dependent) mass 

m« = (4.14) 

which slowly converges to the screening mass defined as lim 2 ^ 00 m(z). At finite density 
this obviously does not work, because C z vanishes at z ~ 1//j. We have chosen to define 
the mass strictly as the coefficient of the exponential falloff, while keeping in mind that the 
correlator is generally more complicated than in the p = case. 

4.2 Comparison with lattice 

At vanishing quark chemical potentials the high temperature mesonic screening masses 
have been measured on the lattice twenty years ago [15,30]. Recent results can be found 
in [20,21]. They have already been discussed more thoroughly in [19]. The masses are found 
to be close to the ideal gas limit, with ir meson lying somewhat below p. The deviation 
from the free theory result in those lattice computations is invariably negative, with masses 
around 5-10% below 2ttT, whereas our NLO analytic calculation gives screening masses 
slightly above the noninteracting result. The results of the current paper behave likewise, 
since the p = limit of our result agrees with our previous paper. 

Lattice measurements of hadron screening masses at finite density have been performed 
only recently, mostly because of the complex action which makes the simulations difficult. 
In [22] the first two coefficients in the Taylor expansion of the mass around p = are 
measured for both isoscalar and isovector chemical potentials. The simulations are done 
at /j = 0, and the result is of course valid only when the chemical potential is small. The 
leading term is just the screening mass at p = 0, and again this lies a little below the free 
theory value, ir (pseudoscalar) below p (vector). The first derivatives vanish at p = 0; for 
ps this follows from a symmetry, and for py it is measured to be consistent with zero. 
This agrees well with our results. 
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The /i-dependence of the screening masses is given by second derivatives measured 
in [22]. They observe that for both ir and p the second order response increases with 
temperature in the isoscalar case, while in the isovector case it approaches zero from below. 
Note that the simulations in [22] were done with two dynamical quark flavors, corresponding 
to dashed lines in Fig. ||. The isoscalar case qualitatively agrees with our result, with the 
mass increasing roughly quadratically near pg = 0. Also in the isovector case the slowly 
decreasing mass in Fig. |5] agrees with the small negative second derivative in [22]. It 
should be noted that at least in our analysis the small value of the second derivative is a 
consequence of the choice Np = 2, giving the dimensionless parameter p in Eq. ( 4.10 ) as 

l + p 2 



P 



(J>v 
vrT' 



(4.15) 



which depends only weakly on py. F° r Np = 3 the derivative vanishes completely, and 
Np = 4 would give a positive second derivative also in the isovector case. 

To get some quantitative comparison, we can also try to fit a quadratic curve to our 
data. At larger p/nT the behavior of our result is clearly something else than simple p 2 , 
so we quite arbitrarily restrict the fitted region to p/ttT < 0.5. Assuming that the first 
derivative vanishes, we look for the coefficients of Re(-Eo) = c i + C 2(p/^T) 2 , which gives 
the derivatives as 



d 2 Re(£ ) 



-0.62, (N F 
0.45, (Np 
0.75, (Np 
0.34, (Np 



0) 
2) 
3) 

3,/x s 



d 2 Re(£ 
dp 2 , 



0) 



-1.42, (JV> 
-0.31, {Np 
0.00, (Np 
-0.41, (Np 



0) 
2) 
3) 

3,yUs 



for isoscalar and isovector, respectively. In terms of physical parameters 

d 2 m _ g 2 C F d 2 Re(E ) 



d,u 2 2tt 3 T dp 



r,2 



-0) 
(4.16) 



(4.17) 



which depends on g^. Using g\jT ~ 2.2 as we did earlier, this gives second derivatives 
of order Td 2 m/dp 2 ~ ±0.02. For isoscalar they are much smaller than in [22], but the 
difference may partly follow from the different definitions of mass. For example, including 
the cosine term like in Eq. ( 4.13j ) would contribute an additional 2T/M ~ 1/ir. 



5. Conclusions and outlook 

This paper has extended our previous computation of mesonic screening masses to sit- 
uations where not only high temperature but also a finite chemical potential is present. 
For flavor non-singlet operators we have computed the full correlator to leading order in 
perturbation theory and derived a next-to-leading order correction to the screening mass. 
We have studied how the mass depends on the chemical potentials of the quarks as well as 
on the number of dynamical flavors. 

Our method is valid for an arbitrary combination of chemical potentials. The numerical 
value of the NLO correction was computed for two distinct cases. We found that when 
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Np > 0, the screening mass increases with jjl if both quarks in the operator have the same 
(isoscalar) fx, whereas for opposite chemical potentials the mass slowly decreases. With 
quenched quarks the mass decreases in both cases. The data near ji = is well described 
by the quadratic curve m(/x) = m(0) + Cfi 2 , but the coefficients we get for the second-order 
response are an order of magnitude smaller than in the recent lattice data. 

The meson correlator looks like a simple exponential only at asymptotic distances. At 
shorter scales the contributions from higher Matsubara modes as well as from factors of 1/z 
modify the behavior, as can be seen from plots of the effective mass. At finite density the 
situation is even more complicated, because the chemical potentials manifest themselves as 
oscillatory terms, which can have a large effect on the short distance behavior. In fact, at 
scales where the correlator is significantly different from zero the screening mass alone does 
not give a sufficient description. It has its virtue in being a simple parameter characteristic 
of screening phenomena, accessible both analytically and on the lattice where there is a 
long history of measuring the screening correlator. However, it is not easy to measure 
the correlator on the lattice at large distances because of the exponential vanishing of 
the correlator and large lattices needed. On the other hand, to extract the mass from 
short-distance data we need to know the general behavior of the correlator. 

We have not touched the issue of flavor singlets in this paper. At large distance their 
correlators are dominated by couplings to different glueball states, whose masses are lower 
than 2irT. We have worked out the couplings in [19], and the glueball masses have been 
measured in [13] also for finite chemical potentials. One notable change in finite density 
is that the Tt^q operator corresponding to V^oV' now couples to TcAq, which has a much 
lower mass. 

Even the finite chemical potential does not change the fact that our NLO correction to 
the screening mass is strictly positive, unlike all the lattice measurements. At asymptot- 
ically high temperature the perturbation theory should give correct results, so we expect 
the screening mass to go above 2ttT as temperature is increased. On the other hand, imple- 
menting our theory on the lattice could show if the higher order and nonperturbative terms 
bring the result below the free theory value at temperatures close to the phase transition. 
The computed correction is small even when g is large, but the next linear 0(g 3 ) term in 
the potential could easily be of the same magnitude. Our theory could also be used to ex- 
tract other features of the correlators besides masses. This would require matching also the 
fields to g 2 order and probably solving the Schrodinger equation in some approximation. 
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A. Correlators in free theory 
A.l Isoscalar chemical potential 

For isoscalar chemical potential /U u = = or for flavor singlets the correlator C z [O a , O b ] 
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in Eq. ( 2,10 ) can be calculated explicitly in free theory. The details depend somewhat 
on the spin structure of the operator, we give here the computation for the scalar (and 
pseudoscalar, up to a sign): 



C z [S a ,S b } = I dr Jd 2 x ± (TP(T,x ± ,z)F a iP(T,x ± ,z)^(0,O)F b iP(0,Q)) 

i r l / T r 

= --8 ab J dr yd 2 x ± (^(r,x ± ,z)V;(0,0))^(0,0)^(r,x ± ,z)) 



i/T r- f TrM f f 

" x * ' £?_ i{u}{ n -u} J n )T+i(p ± -q ± )-x.±+i(p3-q 3 )z 



ipq p 2 q 2 

70 d^3 (Wn -iflf +Pi+ P3Q3 



= 2N c T5 ab S" / / — ^rt-W TPl_^PM3 e i( PS -q 3 )z_ 

V Jp J-°° 27T [(^ - W) 2 + v\ + pj] [{uL - w) 2 + Pi + 9§] 

The chemical potential and Matsubara frequencies only appear in terms like (ojn — i/x) 2 , so 
we can restrict the summation to positive modes by adding the term with \x — ► —fi. The 
poles in p% and (73 can be found at 



P3, <?3 = ±iy (wn - i^) 2 +P 2 L = ±w n (P±), (A.l) 

where the square root is well defined off the negative real axis, and we choose the branch 
with positive real part. Integration over p% and r/3 then gives 



n>0 

12 



2vr 

n>0 

= jab d j^cT 1 e 2 ifiz + / _^ _ x = ^afe d N cT COS 2^Z 

dz 87rz 2sinh27rTz dz &irz sinh27rT,z 

NT/ 1 \ 

= <5 afe \ 2vrT coth 2vrTz cos 2fiz + 2/i sin + - cos 2^z . (A.2) 

87rz sinh 2irTz \ z ) 

A.2 Arbitrary chemical potential 

When no restrictions are imposed on the chemical potentials appearing in the correlator of 
Eq. ( p. 10 ), even the free theory result becomes very complicated. Proceeding in the same 
way as we did above leads to 



C,{S°,S b ] = F$Fy\—^— (e- A ""' T - 1 



m,n 



d 1 



A/iy + i2-7rT(m — n) dz z 
- (A//ij + i2-KT(m - n)) Hdz'^e-^^-^+y/^-wW* 

J z z 
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where A//;j = //i — fly We will also use the shorthand notation p,^ = m + //j below. 

The dominant contribution to the Matsubara sum comes from the lightest modes 
u> = ±7rT, with the remaining sum suppressed by exp(— 2-kTz). Summing only over the 
first modes is then a good approximation to the whole expression. We drop those subleading 
terms and get 

C Z [S\ S b ] = F^F^- (e-^' T - l) |- 



\ r d i 



+ 



■(A/Ujj cos A^ijz + 27rTsin A^y s z 



COS /Uij z 

-2ttT2; 



n -27r7V 



(A/iy)2 + (27TT)2 

— / — - [A/Xij cos /iijz' + A/Ujj cos A/Xijz' — 2-7rT sin A^yz'] e 

11 

(A/ijj cos Afiijz + 2irT sin A/iyz 



+ sin A/xyz + 



-2ttTz 



(A.3) 



z(A / u ij )2 + (2^T)2 

+ A^ijRe[Ei(-27rTz + i/Zyz) + Ei(-27rTz + iA^jz)] - 27rTImEi(-27rTz + iAmjz)\. 

To better understand the long distance behavior of this correlator we can approximate the 
exponential integrals by 



Ei(x) = — ( 1 + - + 0(x 



(A.4) 



which leads to our final expression 



a rpb Jv c 



C z [S a ,S b ] =F?F! 



N r T 2 



2vr 



1 - e 



+ 
+ 



z(A Wj )2 + (2vrT)2 
1 1 

+ 



B. Computation of the potential 



-A Wj /T^j -^ Tz (fix sin^z + (2vrT + cos/Zyz) 

(A/Xij cos A/Uijz + 27rTsin A^ijz) 
cos /Zijz — jLfcij sin /2ijZ 

[2/Zij2-7rTsin / UijZ — ((2ttT) 2 — p,Z) cos/iyz 



. (A.5) 



In this Appendix we show how to compute the static quark-antiquark potential in the 
effective theory. To make the notation lighter, we will assume below that we are comput- 
ing with modes oj n > 0. We need to introduce a point-splitting and study the modified 
correlator 

C(r,z) = J d d - 1 R(4(R,z) Xj (R,z) X ](r/2,0) ( /)i(-r/2,0)). (B.l) 
Going to the momentum space and performing the trivial integration over R, this becomes 
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The leading order diagram of Fig. 0(a) consists of two free quark propagators, 



C<°>(r,. 



-2N r 



2N c 6{z) 



d d p dq 3 



d 2 p 



■ cxp 



2p i J 

2 2 

M; + Mj + ^- + ^- 



= i(53-p 3 )z-ip-r 



Z — zp • r 



0( z )AT c M exp 



(B.3) 



where j?oij is given by Eq. (|4.5| ) and My = Mi + Mj. The high T limit of this expres- 
sion together with its complex conjugate correctly reproduces the approximate 4d result 
Eq. fl2,12|) . In the isoscalar case we have checked that the exact free theory result Eq. ( 2.11| ) 
is recovered if we add the relativistic kinematic correction operator — V^_/8M 3 . 

The NLO diagrams in Fig. 0(b,c) are computed the same way, the only complication 
being the more involved pole structure in p%, q% integrations. For A$ exchange the diagrams 
give 



0(b) = slC F 



0(c) = sICf 



d d k 



(2n) d k 2 + 
d d k 



1 ^(l-e ih A C<®(r,z) 
- mt kn V / 



E ^3 
1 



1 



Jk^z 



(2ir) d k 2 + m| k 2 
and similarly for ^3. The full NLO correction is then 

d d k 1 1 



-ik^z 



C(o) ( 



(B.4) 
(B.5) 



c«(r„ 



r,z 
1 1 



(2^) d fc 2 + m 2 k 2 



2 _ e ifc 3Z _ e ~ik-iZ 



2 — e ifc 3^ _ e ~ik 3 z 



1 + e 



kr 



1 



/t 2 + A 2 fc| 

where A is again the infrared regulator in the gluon propagator. 
The leading order term satisfies the equation 



(B.6) 



d z + M« 



1 



2R)i : 



-V 



C<®(t,z) = 2N c 5(z)S 2 (r), 



(B.7) 



whereas the similar equation for the first order term can be formulated as 



8 Z + Mi 



Vl ^ ^ r, z) = - 9 lC^K ( ^, ™, A M, M, }rUn j C (0) (r; z)) 



2p 



2w n w n w n u; n u n uv n 



(B. 



where /C is some dimensionless kernel. The right-hand side is already multiplied by gL, so 
we need to keep only terms of order unity inside the kernel to get the potential to order 



g 2 . According to the power counting rules of Sec. ^2, we can then expand in the first two 
parameters, taking only the leading term into account. This gives 

V 2 



d z + Mi 



2poi 



r. : ! =a -g E C F u n K, ( 0,0, , — , — -, — l ,ru n 

LU n U) n U) n LO n 



(B.! 



- 20 - 



where 



U! n U! n LO n LO n 



lim d z 

z— >oo 



lim 

2— »00 



d d & 1 
(2vr) d A| 



d A; 2sinfc3Z 
(2vr) d A 3 



k 2 + 



nil 



l + e 



ikr 



1 



/fc 2 + A 2 



„ikr 



fc 2 + mi. 



(2ir) d 
d 2 ' 2e k 
(2vr) 2 - 2e 



2vr5(fc 3 



1 



k 2 + A 2 



fc 2 + 



1 



77? 



l + e' 



l + e 1 



k r 



kr 



1 



fc 2 + A 2 



ikr 



E 
ikr 



1 



k 2 + X 2 



ikr 



k 2 + 



l + e 



ikr 



The signs here co-operate just right to give both ultraviolet (e 
finite expression. The integrals needed here are 



(B.10) 

0) and infrared (A — > 0) 



/ 



q2-2e 



k 



1 



-2c 



(2vr) 2 - 2e A 2 + A 2 
d 2 k e ikr 



I + 21n^ + 0(e) 
47r V e A 



(2^/,* + A* -i K0(Ar)= 27(- ln Y-^ + O(Ar 



(B.ll) 
(B.12) 



rriEr 

In— h 7B - Ko(mEr) 



so we have 

y(r) = 9e C F /<C 0,0, , — — ±,rw n = — 

V w n W n W n U n / 27T 

(B.13) 

and the full correlator C(r,z) = C^°'(r,z) + C^\r, z) satisfies Eq. (|4.4j ) as required. Note 
that the expansion in 1/z above corresponds to the static potential t —* oo in the corre- 
sponding 2+1 -dimensional theory. Similar methods have been recently used to compute 
the real-time static potential for a heavy quark-antiquark pair in 4d QCD [31]. The po- 
tential is infrared finite to this order, so we do not need to consider the contribution of the 
non-perturbative, long distance massless gauge fields. 
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